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ABSTRACT 



We study the statistical properties of the wide binary population in the Galaxy field with projected separations larger than 200 AU by 
constructing the stellar angular two-point correlation function (2PCF) from a homogeneous sample of nearly 670 000 main sequence 
stars. The selected stars lie within a rectangular region around the Northern Galactic Pole and have apparent r-band magnitudes 
between 15 and 20.5 mag and spectral classes later than G5 (g — r > 0.5 mag). The data were taken from the Sixth Data Release of the 
Sloan Digital Sky Survey. We model the 2PCF by means of the Wasserman- Weinberg technique including several assumptions on the 
distribution of the binaries' orbital parameters, luminosity function, and density distribution in the Galaxy. In particular, we assume 
that the semi-major axis distribution is described by a single powerlaw. The free model parameters - the local wide binary number 
density «wb and the power-law index A of the semi-major axis distribution - are inferred simultaneously by least-square fitting. We 
find the separation distribution to follow Opik's law (A = 1) up to the Galactic tidal limit, without any break and a local density of 
5 wide binaries per 1 000 pc 3 with both components having spectral type later than G5. This implies that about 10% of all stars in 
the solar neighbourhood are members of such a late-type wide binary system. With a relative statistical (2cr) error of about 10%, 
our findings are in general agreement with previous studies of wide binaries. The data suggest that about 800 very wide pairs with 
projected separations larger than 0.1 pc exist in our sample, whereas none are found beyond 0.8 pc. Modern large-scale surveys make 
the 2PCF method a viable tool for studying wide binary stars and a true complement to common proper motion studies. The method is, 
however, seriously limited by the noise from optical pairs and the (over)simplifying assumptions made to model the selection effects 
and to interpret the measured clustering signal. 

Key words. Binaries: general - Binaries: visual - Galaxy: stellar content - Stars: statistics 



1. Introduction 

Binary stars have traditionally been central for astronomy, espe- 
cially close binary systems (with typical orbital periods of days 
to years), because they are genuine laboratories of stellar evo- 
lution and its exotic remnants, making them a cornerstone for 
determining masses, distances, and many other fundamental as- 
trophysical parameters. On the other hand, wide binary systems 
(orbital periods of many thousands to millions of years) are of 
particular interest, too. 

Because they are only weakly bound by gravitation, wide 
binaries are prone to tidal disruption by passing massive ob- 
jects, such as massive stars, molecular clouds, MACHOs (mas- 
sive compact halo objects), or dark matter (DM) substruc- 
ture. The shape of their separation distribution, in particular to- 
wards the most extreme, widest separations, should therefore 
allow const raints on the mass and frequency of the disruptive 
perturbers dRetterer & King 1982I). as wel l as to estimate the 
age of a population (Poveda & Allen 2004). Wide binary -based 
MACHO constraints have been placed (but subsequently criti- 
cised) by several auth ors (Bahcal l et al.lll985at [W einberg et al. 
19 O|Weinberg|[T990llWasserman & Weinberg 99 UlYoo et all 
20041 louinn et all2009h . ~ 



In a different context, iHernandez &~Leel |2008) propose us- 
ing the tightening of wide binaries in dwarf spheroidal galax- 
ies through dynamical friction as a test for DM. Moreover, bi- 
naries with separat ions over 0.1 pc, which ar e known to exist 
in the Galaxy (e.g. Lepine & Bongiorno 2007, also this paper), 



are in the 'weak-accelaration' regime where, in principle, one 
could test f or possible deviations from Newtonian gravity, such 
as MOND (iMilgrom & Bekenstein|[l987t IClose et al.lll990h . 

Wide binaries are not only probes of dynamical evo- 
lution and Galactic structure, but they a lso p rovide im- 
portant clues to star formation (Larson] 1200 ll) . The ex- 
act outcome of the binary population in a star-formation 
event is an exceedin gly complex and still unsolved problem 
dGoodwin et al.ll2007l and references therein). The formation of 
extremely wide binaries is particularly difficult to understand 
(lAllen et alj2007tlParker et al]2009t iKouwenhoven et al.l2010l) . 
A good knowledge of the wide binary frequency and separa- 
tion distribution is primary for a whole host of problems in 
astrophysics (e.g. lChanamell2007h . 

Unfortunately, the long orbital periods of wide binary sys- 
tems make their identification very difficult in the first place. 
There are basically two different methods for detecting a wide 
binary system: (1) by the number excess of neighbours around 
a given star with respect to a random distribution, or (2) by 
the common proper motion (CPM) of two well-separated stars. 
Although there is no way to observe the orbital motion of a 
wide binary pair, the CPM is nevertheless a reliable indica- 
tor of a physical relationsh ip between two individual stars (e.g. 
Lepine & Bongiorno 2007). In the present study we use the an- 
gular two-point correlation function (2PCF) to measure the ex- 
cess of neighbours compared to a random distribution. This 
method has the advantage that larger samples of more distant 
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stars can be used, as only the stars' positions are required. It 
only allows, however, statistical statements on the pairing and 
is limited to relatively small angular separations, because the 
noise due to randomly associated pairs increases rapidly (lin- 
early) with angular separation. 

To date, our knowledge on wide binaries essentially had 
rested on studies using the CPM method. T here are two excel- 
lent re cent CP M studies of wide b i naries by IChaname & Gouldl 
(2004 and ILepine & Bongiornol (120071) . Based on updated 
proper motion (Luyten and Hipparcos) catalogues of stars, 917 
and 521 wide binary systems, respectively, have been identified 
over the entire N orthern sky and with a typical median distance 
of about 100 pc. ILepine & Bongiornol found the separation dis- 



tribution of the pairs to follow Qpikl's (1 19241) law, i.e. frequency 
being proportional to the inverse of separation, out to separations 
of around 3 500 AU. Beyond this characteristic scale, however, 
the separation distribution seems to be falling by a steeper power 
law, without an obvious cut-off. ILepine & Bongiornol find wide 
binaries out to separations of almost 100 000 AU, or 0.4 pc! They 
quote a number of 9.5% for nearby (D < 100 pc) Hipparcos stars 
belonging to a wide binary system with a separation greater than 
1 000 AU, agai n demonstrating th e ubiquity of the (wide) binary 
phenomenon. Cha name & Gouldl in addition, achieve a distinc- 
tion between wide binaries belonging to the Galactic disc and 
those belonging to the Galactic halo. No significant difference 
in the separation distribution has been found, suggesting that the 
disc and halo binaries were formed under similar conditions, de- 
spite the very different metallicities and ages. 

The angular 2PCF (also called "covariance function") is 
one of the most useful tools for studying the clustering prop- 
erties of galaxies, and it has been widely used to prob e the 
large-scale structure of the Universe (e.g. iPeeblesI I1980L for 
a more recent study o f the angular clustering of galaxies 
see lScranton et alJl2002l and references therein.). The 2PCF can 
also be, and has been, used to measure the clustering proper- 
ties of stars. While on large scales the stars are clearly randomly 
distributed in the Galaxy, as expected in a collisionless system, 
correlations up to the 10 pc scale (or ev en beyond) can be foun d 
in moving groups and hal o streame rs (Doinidis & Beers! ri989h . 
star-formi ng regions (e.g. [Gom ez & Lada 1998), and open star 
clusters (lLopez-Corredoira et alj U998). On very small scales 
(sub-pc, observationally: sub-arcmin), there is a strong signal 
due to (visual) binary stars. 

The 2PCF method of pro bing wide binary stars i n the 
Galactic field was pioneered by Bahc all & Soneiral d 198 lb . who 
studied the distribution of stars down to a limiting magnitude V 
= 16 in a 10 square degree field at the NGP, and found very 
significant clustering at a separation of 0.1 pc. Of the 1 9 bi- 
nary candidates, 6 turned out to be real (lLathametalJl 984). The 
theoretical implications of these observations for the frequency 
and separation distribution of the binari es, and a general method 
for mo delling them, was worked out by Was serman & Weinberg! 
(fl987h . To da te, there are only few follow-up studies of the 
stellar 2PC F (Saarin en & Gilmordl 1989tlGarnavichlll988l.il 99 U 



lGouldetalJll995l) . Given this surprising lack of further work on 
the stellar 2PCF to study wide binaries, and in view of the enor- 
mous progress in deep photometric sky surveys in the past 20 
years that should render the 2PCF method - as a true comple- 
ment to the CPM method - much more rewarding now, we have 
started a project to use t he huge stellar d atabase of the Sloan 
Digital Sky Survey (SDSS lYork et al.l2000l www . sdss . org) for 
a detailed stellar correlation analysis to very faint magnitudes. 

An indepen dent study of wide binaries in the SDSS database 
bv ISesar et al .] (I2008h takes a different approach, based on the 



"Milky Way Tomography" of lJuric etal 1 (12008b . where approxi- 
mate distances are ascribed to all stars by adopting a photometric 
parallax relation. Candidate binaries are selected by the require- 
ment that the difference between two potential components in 
apparent magnitude is within a certain error equal to the differ- 
ence in absolute magnitude. 

Using distance information has the fundamental advantage 
of filtering out the disturbing noise fro m chance projections and 
of avoiding the need for a complex, Wasserman & Weinberg- 
type modelling to calculate integrated, sky-projected quantities. 
Although we plan to include distance information in future work, 
we show here that the angular correlation analysis is in principle 
still a viable approach. Our main results agree with those of the 
CPM studies. We also show, however, what the limitations of the 
method are. 

The paper is organised as follows. In Sect. [2] we describe 
the SDSS data input and define the cleansed sample used for 
the analysis. Section [3] is devoted to the angular correlation 
function apparatus, followed by an extensive description of the 
modellin g of the correl ation function drawing on a modified 
IWasserman & W einberg technique in Sect. [4] Our results for the 
total sample and for a number of subsamples are presented in 
Sect. |5] and are critically discussed and compared to previous 
work in Sect. [6] Concluding remarks are given in Sect. Q 



2. Data 

For our analysis we selected a homogeneous sam- 
ple of stars from the Sixth SDSS Data Release (DR6; 
lAdelman~ McCarth v et al.l 120081) . We took a rectangular area 
around the Northern Galactic Pole (NGP) covering approx- 
imately Q tot ^ 675 square degrees (dec: 22° - 42°, RA: 
165° - 205°; see Fig. [T}. It contains 966 656 primar$\ point-like 
objects - including quasars, asteroids, and possibly some 
misidentified galaxies - selected from the Star viewQ, having 
an apparent PSF magnitude in the r-band between 15 and 
20.5 mag. Following the SDSS recommendation^ we required 
the stars to have clean photometr^ in the g, r, and i band. 
Although the effect of interstellar dust on the measurements is 
wa ek in the direct i on of the NGP, we corrected the data using 
the ISchlegel et alJ (ll998h maps, which can be easily done, as 
the extinction at the position of each object is stored in the 
SDSS database. 

The magnitude limits were chosen to be within the satura- 
tion l imit of the SDSS CCD cameras (r ~ 14 mag: iGunn et al.l 
ll998lFI and the limit whe re the star-galaxy separation becomes 
unreliable (r ~ 21.5 mag: lLupton et alJl200llrl Additionally, we 
avoid close stars being overcorrected by the adopted extinction 
correction, since most stars with r > 15 mag are behind the en- 
tire dust column dJuric et al . 2008). On the other hand, with the 
somewhat conservative faint limit of r — 20.5 we make sure that 
only very few stars have a magnitude in the g or i band beyond 



1 Due to overlaps in the imaging, many objects are observed more 
than once. The best of those observations is called the primary object, 
all the others are called secondary objects. 

2 See cas . sdss . org/ dr6/ en/help/browser/browser . asp 

3 See sdss.org/dr6/algorithms /photometry. html#mag_psf 

4 See sdss.org/dr6/products/catalogs/flags.html 

5 See cas . sdss . org/dr6/en/help/docs/realquery . asp#f lags 

6 See also astro . princeton . edu/PB00K/camera/ camera . htm 

7 See also sdss . org/dr6/products/general/stargalsep . html 
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Fig. 1. Distribution of stars (grey points filling the background) in our total sample. Black asterisks show the positions of bright star 
masks, white circles those of hole masks. The size of the symbols are not to scale. 



the 95 % completeness limit of 22.2 m ag or 21.3 mag, respec- 
tively dAdelman-McCarthv et ai]|2007l Table lfl 

The average seeing of the SPSS imaging data (median PS F 
width) is 1.4" in the r-band (lAdelman-McCarthv et all 12008 fl 
To be on the safe side, we took the minimum angular separation 
to be min = 2". 

2.1. Contaminations 

Matching our sample with the QsoBest tabled resulted in the 
exclusion of 10041 quasar candidates. Most of them (8157) 
have g - r < 0.5 mag and are scattered in a colour-colour di- 
agram around the otherwise narrow stellar locus as shown in 
Fig. |2 Even after removing the quasar candidates, the remaining 
objects in the blue part of our sample show a suspicious scat- 
ter, which is probably caused by further quasars and misidenti- 
fied galaxies. We therefore decided to exclude all objects with 
g - r < 0.5 mag, removing a further 286 227 objects from our 
sample. Furthermore, we reject all moving objects (asteroids) by 
cutting on the DEBLENDED_AS_MOVING flag. This leaves us with 
670388 objects classified as "stars" in the sample. 

The large majority of th e stars observed by the SDSS are 
main sequence (MS) stars. Finlator et al.l d2000h estimate the 
fraction of stars that are not on the MS to be ~ 1 .0%, most of 
them giants and subgiants (~ 90%), but also horizontal branch 
stars (~ 10%). The number of white dwarfs observed by the 
SPSS is negligib le compared to the number of MS stars (e.g. 
lHarris et alJl2006l) . Thus, it seems well-justified to assume that 
all the stars in our sample are on the MS. 

The cut discussed above at g — r = 0.5 mag implies that our 
sample contains only stars from spectral type later than about G5 
(Finl ator et al.ll200U) . In addition, this cut is appropriate for our 
purposes for the following three reasons: 

8 Seealsosdss.org/dr6/ 

9 See also http://www.sdss .org/dr6/ and the DR5 pa- 
per (|Ad^^McCMftyiral]|2003) 

10 See cas . sdss . org/dr6/en/help/docs/ algorithm . asp?key= 
qsocat 



i) Because they are very young, the bluest MS stars are 
mostly members of loose associations, so their clustering proper- 
ties still represent the peculiarities of their birth places. Bei ng an 
interesting subject to study (e.g. lKobulnickv & Frverll2007h . ex- 
cluding them assures that our clustering signal is predominantly 
due to wide binaries in the field that have lost their memory of 
their birth places. 

ii) As the MS becomes more sparsely populated towards the 
blue end, and a significant fraction could be made up by metal- 
poor halo giants, the assumption that all stars in our sample are 
on the MS might not be valid for the bluest stars. 

Hi) For magnitudes My > 4.5 mag, the shape of the halo 
luminosity fun ction agrees well w ith that of the disc luminosity 
function (e.g. iBahcall et al]|1985bi their Fig. 2). Thus, the cut 
on g - r = 0.5 mag, which corresponds to a cut at My — 5.6 
mag, allows us to use the disc luminosity function for the halo 
component. 

2.2. Survey holes and bright stars 

The sample chosen contains some regions where, for vari- 
ous reasons, no object could be observed. These regions are 
maskecQ as hole in the SPSS database and can therefore be 
easily identified. Holes in the sample affect our analysis in two 
ways: first, they diminish the total area of the sample. Second, 
they introduce an edge effect, as stars near a hole show a lack 
of neighbouring stars. The former effect can be easily corrected 
for in an approximate way as the SPSS provides the radius flj^ 
of its bounding circle for every mask i. The residual area of the 
sample is then 

Q-ft tot -£7r(0®) (1) 

i=\ 

where Nm denotes the number of hole masks. We discuss the 
correction for edge effects in Sect. l3.2l in the context of the cor- 
relation function. 

11 See sdss.org/dr6/algorithms/masks.html 
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Fig. 2. Colour-colour diagram of all point-like objects within 
our coordinate's and apparent magnitude's limits (grey points). 
Black points are SDSS quasar candidates. The vertical dashed 
line indicates the cut at g - r — 0.5 mag (see text). Approximate 
spectral classes are indicated. 



Restricting ourselves to the masks defined in the r-band, we 
find A^sh = 9 regions masked as "hole" in our sample (survey 
holes: SH) with an average radius of Ssh - 7.76 arcmin. These 
survey holes diminish the total area of the sample by approxi- 
mately 0.46 square degrees, constituting a marginal correction 
that could be safely ignored. 

Very bright (saturated) stars cause similar problems: a very 
bright object may appear like a hole in our sample, when the 
underlying fainter stars blended with this object cannot be re- 
vealed. Within our sample there are A^bs = 1 790 bright star 
(BS) masks in the r-band. We exclude all objects inside such a 
mask, removing 545 stars from our sample. The bounding cir- 
cles of bright star masks have an average radius of #bs — 1 -95 
arcmin, resulting in a further diminishing of the total area of the 
sample of circa 5.96 square degrees. By correcting the sample's 
total area for both hole and bright star masks, we get a residual 
area O ^ 675 - 0.46 - 5.96 deg 2 = 668.58 deg 2 . 



2.3. Final sample 

Our final sample contains N b s = 669 843 MS stars with clean 
photometry, a spectral type later than about G5 (g-r > 0.5 mag), 
and an apparent magnitude in the range 15 < r < 20.5 mag. 
The stars are distributed over a solid angle Q 668.58 square 
degrees (after correcting for bright star masks and minor survey 
holes), corresponding to a mean surface density of n — No^/Q. - 
1 000 stars per square degree. Throughout this study, we use the 
terms "total sample" and "final sample" synonymously. 



3. Stellar correlation function 

3.1. Estimation of the correlation function 

The angular two-point awfo-correlation function (2PCF) w(9) is 
defined via the joint probability dP of finding objects in both the 
solid angles dQi and d£^2 



dP = « 2 (1 + w(0))dQ,dQ 2 , 



(2) 



where n is the mean density of object s in the sky a nd 9 is the 
separation between the two areas fe.g. lPeeblesl [T980. §45). The 
number of distinct stellar pairs, F(9)d9, with an angular separa- 
tion between 9 and 9 + d9, observed in a region of angular size 
£2, can then be written as 



F(9)d9 = 



N obs (N obs ■ 



1) dO(0) 
— (1 + w(9)) ■ 



n 



(3) 



where A^bs is the number of stars observed in O and d£l(9) is 
approximately equal to the solid angle of a ring with (middle) 
radius 9 and width d9 

d£2(60 - 2n6de . (4) 
Solving for w in Eq.[3]yields a simple estimator for the 2PCF: 



P(9) 



(5) 



where P(9) is the number of pairs expected from a random sam- 
ple with A^obs points (w = in Eq.0: 



P(9)d8^ 7TNoUN ^- 1) 9d9. 



(6) 



The 2PCF estimate w(ff) is a measure for the excess of observed 
pairs separated by an angle 9 with respect to a randomly dis- 
tributed sample. 

Another statistical measure equivalent to the 2PCF, which 
we use for visualising the data, is the cumulative difference dis- 
tribution (CDD) y{9), giving the cumulative number of pairs in 
excess of a random distribution: 



y{9)= {F{ff)-P{ff))dff 



w(9')P{ff)d8' 



(7) 



with flmin = 2" (see above). The CDD is completely equivalent to 
the 2PCF, but is - as for all cumulative distributions - more sen- 
sitive to statistical trends caused over a wider range of angular 
separations (i.e. over several bins), at the expense of a strong cor- 
relation of its values at different separations. The use of the sim- 
ple estimator in Eq.[5]is appropriate as long as b oundary e ffects 
due to the finite sample size are negligible (iBerns tein 1994). We 
discuss boundary effect in the next section. The number of pairs 
observed F{9) is determined very efficiently inside the SDSS 
database using the precalculated Neighbors table which contains 
all the objects within # raax = 30" of any given object. 

3.2. Boundary effects 

Stars close to the sample's boundary have a somewhat trun- 
cated ring dO(#) and consequently show a lack of neighbours. 
However, as in the present study the probed angular scale is 
small compared to the sample's size, it turns out that these edge 
effects have a negligible impact on our results: we estimate the 
relative error introduced in dQ.(9) by omitting the edge correc- 
tion to be less than ~ 0.04% for any given 9 < 30". 
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The stars near the boundary of a hole or a bright star mask 
(in the following simply "hole") show a lack of neighbours, too. 
Therefore, only a fraction F™(0) of all pairs separated by an 
angular distance 9 has been observed, and we need to correct 
the observed number of pairs F{9). The "true" number of pairs, 
corrected for edge effects, then reads as 



F C on(9) - 



F(9) 
F™(9) 



(8) 



In calculating the fraction F™ we proceeded in the same way 
as lLopez-Corredoira et alJ (fT99 8). The calculation is outlined in 
Appendix [A] We find that the relative error in F for 9 = 30" 
amounts to 



A F = 



= 1 - F„ a 0.7% 



(9) 



For smaller angular separations, the effect is even less. 

Taking into consideration these edge effects, we rewrite the 
2PCF estimate 



w C0lT (9) = 



FcoAO) 



P{9) 
and the CDD 

7coir(0) : 



- 1 



T (6O^(0')d0' . 



(10) 



(ID 



3.3. Uncertainty of the correlation function estimate 

As the values of vv corr at different separations are not inde- 
pendent, Poisson errors may underestimate the uncertainty in 
Wcorr, especially when t he angular scale s under consideration are 
large (Hamiltoi3 [l993t lBernsteinlll994h . However, as we see in 
Sect. [5] the clustering of wide binary stars occurs on small angu- 
lar scales (9 < 15") compared to the size of our sample (20 deg 
x 40 deg). We therefore expect that using Poisson errors only 
underestimates the true errors by a small amount in our case. 
Thus, we adopt Poissonian errors on F(9): 



8F{0) = y/F(ff) and 6F cm (9) = 



F C on(9) 



(12) 



Using Eq. \TU\ and Gauss' error propagation formula, we may 
write the uncertainty of the CF estimate w con as 



6w c , 



SF a 



vvv, 



+ 1 



ptotp 
H 



(13) 



where we omitted the dependencies on 9 for the sake of brevity. 

The uncertainty in the CDD y is easily obtained in the same 
way: 



<5fcon-(#) 



6Fl on {9'W ■ 



(14) 



3.4. Testing the procedure for a random sample 

To test the validity of the procedure to estimate the 2PCF de- 
scribed in the previous sections, we generated a random sample 
having the same number of "stars" and distributed over the same 
area as the stars in our final sample. In addition, we made certain 



0.05 



D. 

o 




400 



200 



CD 

E 

o 



-0.05 



200 



-400 1 



u 



600 



6 8 10 12 14 16 18 20 22 24 26 28 30 
Projected angular separation in arcseconds 

Fig. 3. 2PCF as inferred from a random sample (solid circles, 
left ordinate) and the corresponding CDD (open circles, right 
ordinate). Poisson errors are indicated as vertical lines. 



that the random sample contains the same number of holes with 
appropriate radii. 

For the analysis of the random sample we wrote a dedicated 
computer program that lays a grid with a mesh size of 1 arcmin 
over the sample before counting the pairs. This grid acts as a 
distance filter and avoids calculating the distances between all 
possible pairs: only pairs with both their components within a 
cell or with them in two adjacent cells, respectively, were taken 
into account. 

If our procedure is correct we would expect a zero signal in 
both the 2PCF and the CDD when analysing a sample of ran- 
domly distributed stars. Figure[3]shows the result of the analysis 
of our random sample. There is indeed no clustering signal ev- 
ident, neither in the 2PCF nor in the CDD. The results are con- 
sistent throughout with a zero signal out to the separation limit 
of 30"(compare also with Fig.|5]l. This shows that our procedure 
for estimating the 2PCF is reliable. 

4. The model 

Our approach to mo delling the angular 2PCF i s base d on a tech- 
nique d eveloped by IWasserman & Weinberg! (1 19871) . hereafter 
IWW87L As we describe in more detail in the following sections, 
this technique makes some simple assumptions on the basic sta- 
tistical properties of wide binaries, and projects these theoretical 
distributions on the observational plane using the selection cri- 
teria of a given (binary) star catalogue and the geometry of the 
Milky Way galaxy (see also Weinb erg!! 19881) . 

Their long periods make it virtually impossible to distinguish 
wide binary stars from mere chance projections (optical pairs) by 
their orbital motion. Therefore, we do not attempt to identify in- 
dividual wide binaries, but we look for a statistical signal stem- 
ming from physical wide pairs in the sample, solely exploiting 
precise stellar posit ion measurements prov ided by the SDSS. 

The original IWasserman & Weinberg! (WW) technique cal- 
culates the projected separation distribution to compare it with 
a sample of wide binaries of known distance and angular sep- 
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aration. In the present study we are only dealing nwith angular 
separations. The calculation of wide binary counts as a function 
of angular separation alone re quires an appropia te modification 
of the WW-technique (see also iGarnavichl fl 99 1 ). which we dis- 
cuss in Sect. 14.31 



4.1. Wasserman-Weinberg technique 

IWW87I developed" a versatile technique for comparing wide bi- 
nary observations with theoretical semimajor axis distributions". 
According to WW87 we may write the number of observed bi- 
naries ifr(s)ds with projected physical separations between s and 
s + ds in a given catalogue of stars as 



ifr(s) ds = n WB Q(s)V(s) ds , 



(15) 



where «wb is the total number density of wide binaries in the 
solar neighbourhood^, Q(s) is the reduced distribution of pro- 
jected separations, and V(s) the "effective volume". 

The total number density «wb is one of the two free para- 
meters in the model that will be inferred by fitting the model to 
the observational data. The reduced separation distribution Q(s) 
contains the physical properties from the wide binaries (semi- 
major axis distribution, distribution of eccentricities, orientation 
of the orbital planes) projected onto the observational plane, 
whereas the effective volume V(s) takes into account the charac- 
teristics of the sample under consideration (covered area, range 
of angular separation examined, magnitude limits), as well as 
the stellar density distribution in the Galaxy and the luminosity 
function. 

As we show in Sect. 14.31 this neat formal splitting in physi- 
cal properties and selection effects will not be possible anymore 
after the required modification on the WW-technique mentioned 
above. At this point, we discuss the two parts - Q(s) and V(s) - 
more in detail. 

4.1.1. Reduced separation distribution 

The reduced separation distribution Q(s) is given by the projec- 
tion of the reduced (present-day) semi-major axis distribution 
q(d) against the sky. Gravitational perturbations due to stars, 
giant molecular clouds, and (hypothetical) DM particles cause 
the semi-major axes of (disc) wide binaries to evolve. Little is 
known about the initial semi-major axis distribution, and usu- 
ally a single po wer-law is assumed, because of its simplicity 
(Weinberg 1988) b ut also because of theoretical considerations 
( Valtonen 1997 0Poveda et al.l2007l and references therein). The 
evolution of the semi-major axis distribution of disc wide bi- 
naries has been modelled by Weinberg, Shapiro, & Wasse rmanl 
( 1987) in terms of the Fokker-Plack equation. Their numerical 
simulations suggest that the semi-major axis distribution evolves 
in a self-similar way for reasonable choices of the initial power- 
law index and the wide binary birth rate function. We decided, 
therefore, to model the present-epoch semi-major axis distribu- 
tion of wide binary stars q(a) by a single powerlaw 



(16) 



9(a) pc 



with the power-law index A, the second free parameter of the 
model. The range in a where Eq. [I6]is a valid description of the 
observed semi-major axis distribution at the same time provides 
the definition of what we consider as a "wide binary star". We 



specify this range by an upper and a lower limit, a m ; n and a max , 
respectively. There is observational evidence that the distribu- 
tio n of the se mi-major a xis of bin ary stars has a break at 0.001 
pc dAbtlll983h . Following lwW87l and lWeinbergl(ll988h . we thus 
take the lower limit, providing the division between the close 
and wide binary populations, to be a m ; n = 0.001 pc (« 200 AU). 
On the other hand, the Galactic tides provide a natural maxi- 
mum semi-major axis aj beyond which no bound orbits can ex- 
ist. Details on the calculation of aj are given in AppendixlBl We 
find it to be of the order of 1 pc. 

We normalise q(a) - hence the name ''reduced semi-major 
axis distribution" - such that 



f 



q(a) da — 1 



giving 



q(a) = cA — \ pc 



where 



ca = i 



(l-^^-fl^]" 1 , for ,1*1 



(17) 



(18) 



(19) 



[In fl T - In Omin] 



for A = 1 



with fl min and aj in pc. 

In a more general treatment, q(a) might also depend on the 
luminosity classes of the bin ary star members, as well as on 
their magnitudes. Like lWW87l we restrict our analysis to models 
where q(a) only depends on the semi-major axis a. 

It is somewhat disputed whether the semi-major axis distri- 
bution has a break at larger a attributed to the disruptive effects 
of the environment on the wid e st bin ary stars. I n their extensive 
work lWasserman & Weinberg! (1 19911) (see also IWeinbergl [l990) 
conclud that "although the data suggest a break in the physi- 
cal distribution of wide binary separations, they do not require a 
break with overwhelmi ng statistical significance", whereas the 
more recent study by Lepine & Bongiornd (120071) shows evi- 
dence for a break at s ~ 3 500 AU (statistically, we have (s) =s a; 
see Eq. l22l . In this context, the question arises whether the as- 
sumption that the data is described by a single powerlaw (Eq.[T6ll 
up to the tidal limit aj can be rejected with confidence. We ad- 
dress this question in Sect. [5] 

Assuming that "the binaries' orbital planes are randomly ori- 
ente d and th at their eccentricities e are distributed un iformly in 
e 2 " dWW87l) we can project q(a) against the sky by (Weinberg 
119881) 



r°° da 

I —q(a)T(s/a) , 

Js/2 



(20) 



where T(x) takes the eccentricity and angular avera ging into 
account for the pairs, whic h may be written as (IWW87L 
Weinberg & Wasserman 1988) 



Ax r^-* [ 

r( X ) = -j o du y 



(u 2 + x)(2 — x — u 2 ) 
u 2 + 2x 



(21) 



With the above assumptions on the orientations of the orbital 
planes and the distribution of eccentricities it can be shown that 
the average observed projected separati on (s) for a give n semi- 
major axis a is basically equal to a (e.g. Yoo et al. 2004]): 



We use the term "solar neighbourhood" for local quantities. 



5n 

< 5 > = — a = 0.98a . 
16 



(22) 
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10 4 10 3 10 2 10 1 10° io 1 



Projected separation s in pc 

Fig. 4. Effective volumes calculated using the parameters that 
correspond to our final sample (Table [2j for the three Galactic 
structure parameter sets as discussed in 34.21 Long dashed line: 
set 1; solid line: set 2; short dashed line: set 3. 



(IWW87h . and assuming further that the stellar density distribu- 
tion p is independent of the stars' absolute luminosities, we may 
write the effective volume as 

where we have the standard relations 

M max (D) = m max -51og 10 (D/10pc) (27) 
M min (D) = m min -51og I0 (D/10pc), (28) 

and p(D) is normalised such that p(0) = 1 . 

In Fig.|4]the effective volumes for the three Galactic structure 
parameter sets (described in the next section) are plotted using 
the parameters that correspond to our final sample (see Table[2]i. 
From the effective volume V(s), it is evident that our study is 
most sensitive to projected separations in 0.01 pc S= s ^ 0.1 pc, 
whereas our study is absolutely insensitive when s S= 0.001 pc 
or s ^ 1 pc. Referring to Eq. [22] we see that the shape of the 
effective volume also indicates that we are insensitive to semi- 
major axis a ^ 0.001 pc and for a X 1 pc, which nicely fits the 
range in a where we assume that the single power-law model 
holds. But it also shows that our analysis is not too sensitive to 
very wide binary stars with semi-major axes larger than 0. 1 pc. 



From the normalisations of q(a) it also follows that Q(s) is nor- 
malised to unity in the range [(s) m i n , (s)max]- 

Substituting rj - yfsjla in Eq. [20] as well as £ = uj V2 in 
Eq.GJJleads to a convenient form for numerical integration 



(23) 



where 



min(l, V^/2«min) '\j^~ r l 2 



C A ( s ) = 2^ C -±JdW Ul f dj 

. / . i 7* \ A 



min(l, Vs/2i3t) 



f + 27/ 2 



(24) 



The integrals in E q. l24l are evaluated using Gauss quadrature 
dPress et al.lll992h . For s far from a m ; n or aj the distribution of 
the projected separations Q(s) is approximately powerlaw 



— 

\pc 



pc 



(25) 



4.1.2. Effective volume 



As we deal with a magnitude-limited sample, we are plagued 
by selecti on effec ts that must be properly taken into account. 
Following [WW871 we do this by means of the effective volume. 
It allows for the solid angle Q covered by the sample, the angular 
separation range [0 m i n , max ] we are examining, and the apparent 
magnitude limits m m \ n and m max . Furthermore, the effective vol- 
ume takes the stellar density distribution p into account, as well 
as the normalised (single star) luminosity function O (both de- 
scribed more in detail in Sect. I4.21 I. 

"Making the plausible assumptions that the intrinsic lumi- 
nosities of the two stars in a wide binary are independent and 
distributed in the same way as the luminosities of field stars" 



4.2. Galactic model 

4.2.1 . Stellar density distribution 

Follo wing earlier work (especially lJuric etaL 2008: iChene"t"aTl 
2001, and references therein), we modelled the stellar density 
distribution of the Milky Way Galaxy by including two exponen- 
tial disc components - a thin and a thick disc - and an elliptical 
halo component whose density profile obeys a powerlaw. The 
contribution of the Galactic bulge is negligible in the direction 
of the NGP and it is therefore ignored in the following. 

The three components are added to yield the total density 
distribution 



p(D) = p(0)p(D) 



with 



P(D) 



p thin (D) + n thick p thick (D) + n halo p halo (D) 



1 +n« 



+ n" 



(29) 



(30) 



which is normalised such that p(0) — 1 . The thick disc and the 
halo component are normalised with respect to the thin disc via 
the normalisation constants « tnick and n hal °, respectively. The 
overall normalisation p(0) is calibrated to produce the observed 
star counts. 

Both the thin and the thick disc popul ations obey a dou- 
ble ex ponential density law of the form (e.g. lBahcall & Soneiral 
1980) 



exp 



h r 



exp 



(31) 



where h r and h- are the scale length and the scale height of the 
disc, z is the object's height above the Galactic midplane, and 
r is its Galactocentric distance in the Galactic plane. We take 
the Sun's distance fro m the centre of th e Galaxy in the Galactic 
plane to be tq = 8 kpc. lJuric et al. I d2008h find that the Sun is lo- 
cated zo — 25 pc above the Galactic midplane dChen et al 
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give zo - 27 pc). It is straightforward to derive the following 
useful relations 



Table 1. Galactic structure parameters 



z = zo + D sin b 
r 2 = r \ + d 2 - 2r d cos i 
d = Dcosb, 



(32) 
(33) 
(34) 



where d is its distance from the Sun in the Galactic plane and 
t and b are its Galactic longitude and latitude, respectively. The 
stellar density is fairly constant within the sample, so we ne- 
glect the direction-dependent density variations. For the sake 
of simplicity, we adopt the coordinates of our sample's centre: 
I = 175.6° and b = 81.6°. (Dividing our sample into subsamples 
and summing over them taking their centres hardly influences 
our results.) 

In line with lJuric et ail d2008l) and lChen et all (120011) we as- 
sume the scale heights to be independent of absolute magnitude 
M. Using capitals for the thick disc's and lower case letters for 
the thin disc's scale height and length, we may write the nor- 
malised number density distribution of the stars in the Galactic 
thin disc as 



~thin/ m r-r z\ zo 

p ( D ) = exp^— --j-exp^- 

and that of the thick disc as 



-thick, ^ 1 r — ro z\ zo 
P (°) = ex P| - Jf J ■ ex P ( Jf 



(35) 



(36) 



where the rightmost factors are for normalisation purposes and 
ensure that p^O) = p MA (0) = 1. 

The density distribution of the stellar halo population is mod- 
elled by a powerlaw with index k. The observational data prefer 
a somewhat oblate halo giving an ellipsoid, flattened in the same 
sense as the Galactic disc, with axes a — b and c = Kg , where k 
controls the ellipticity of the halo (cf . lJuric et al . 2008) 



~halo (D) 



r 2 + 



(if 



(37) 



which, of course, also satisfies p halo (0) = 1 . 

We compare three different sets of structure param eters: 

Set 1: The measured values from lJuric et a l. (2008) (see their 
Table 10). These values are best-fit parameters as directly mea- 
sured from the apparent number density distribution maps using 
their "bright" photometric parallax relation (their Eq. 2). They 
are not corrected for bias caused by, e.g., unresolved stellar mul- 
tiplicity, hence the term "apparent". 

Set 2: The bias-corrected values from lJuric et al ] (120081) 
(see again their Table 10). These values were corrected for 
unrecognised stellar multiplicity, Malmquist bias, and system- 
atic distance determination errors by mea ns of Monte Carlo- 
generated mock cata logues (lJuric et al.ll2008l §4.3.). Following 
iReid & Gizisl (1 19971) . the fraction of "stars" in the local stellar 
population that in fact are unresolved binaries is taken to be 35%. 
However, the halo component was not included in the mock cata- 
logues, and its structure parameters were therefore not corrected, 
but the measured values are used instead. This set of parameters 
will be referred to as our standard set. 

Set 3: As a third independent set of Galactic st ructure pa- 
ramete rs, we refer to the somewhat earlier work of IChen et al.l 
J2001. It is also based on observations obtained with the SDSS, 
but Chen et all (12001 1) infer the density distribution of the stars 



structure 


set 1 


set 2 


set 3 


parameter 








zo [pc] 


25 


25 


27 


h- [pc] 


245 


300 


330 


K [PC] 


2 159 


2 600 


2250 


R thick 


0.13 


0.12 


0.0975 


H z [pc] 


743 


900 


665 


H, [pc] 


3 261 


3 600 


3 500 


n halo 


0.0051 


0.0051 


0.00125 


K 


0.64 


0.64 


0.55 


k 


2.77 


2.77 


2.5 



by inverting the fun damental equation of stellar statistics (e.g. 
iKarttunen et al.Hl 9961) . 

These three set s of Galact i c struc ture parameters are sum- 
marised in Table [Tl lJuric et all (|2008) quote unrecognised mul- 
tiplicity as one of the dominant sources of error in the distance 
determination of the stars. (Only the uncertainties in absolute 
calibration of the photometric parallax relation might be even 
more important, but little can be done to incr ease its accuracy at 
the moment; see, however, Se sar et al. I d2008h .) It is therefore not 
surprising that the scale heights and lengths of the bias-corrected 
parameter set are larger than those of the measured parameters 
as the misidentification of a binary star as a single star results in 
an underestim a tion o f its distance. Regarding the values derived 
bv IChen et al ] (120011) . we find the largest difference in the nor- 
malisation of the halo component th at is about four times smaller 
than that given bv lJuric et all (120081) . whereas the other parame- 
ters are not too dissimilar. 



4.2.2. Stellar luminosity function 

For our study we refer to the l uminosity function (LF) in- 
ferred by JahreiB & Wielenl £1997) using Hipparcos parallaxes, 
which gives reliable values for a wide magnitude range (— 1 < 
My < 19). (The f a int en d of the LF is somewhat uncertain and 
iJahreifi & Wielenl (1 19971) give only a lower limit in the range 
20 < My < 23. We take that lower limit at the faint end to 
be the true value of the LF in that magnitude range.) 

We need to transform the IJahreifi & Wielenl (119971) LF 
from visual (V-band) into r-band magnitudes, i.e. from <E>(My) 
to $}(M r ). We perform this transformation in an unsophisti- 
cated way by combining the photometric parallax relation from 
lLaird et al.l (Il988l) . which was calibrated using the Hyades and 



is linear in (B - V) (their Eq. la) 
M v = 5.64(B- V) + 1.11 mag 



(38) 



with a transfor mation equation that is obtained by subtracting 
Eq. 1 from 3 in lBilir et all d2005l) 



V - r = 0.491(B - V) - 0.144 mag . 



(39) 



The linearity of the photometric parallax relation removes any 
inversion problem, since it assures - in a mathematical sense - 
that a colour index {B - V) exists for every My. For a given 
distance we then have 



M r = 0.913M V + 0.0476 mag 
and the transformed LF is simply 
dM v 



<E>(M,.) = <t>(M v ) 



dM,. 



(40) 



(41) 
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Using linear interpolation (the brightest bin at M r = — 1 mag is 
an extrapolation), we finally have rebinned <S)(M r ) such that the 
bins are centred on integer values of absolute magnitude M r and 
have a width of AM r = 1 mag. 

We normalise the LF by integrating over the absolute mag- 
nitude range corresponding to our total sample 



Table 3. Best-fit values: total sample 



d>(M r ) 



0(M r ) <D(M r ) 



poo 



(42) 



with n' t a 0.094 pc -3 being the total number density of stars in 
the solar neighbourhood having g - r > 0.5 mag. Assuming that 
all stars are on the MS , we find, u s ing th e ("bright") photometric 
parallax relation from lJuric et alj d2008l) . that the cut at g - r = 
0.5 mag corresponds to M, 5.2 mag . (We first transformed it 
from g — r into r - i using Eq. 4 from lJuric et al.h 

4.3. Modification of the Wasserman-Weinberg technique 

To derive a theoretical angular 2PCF we need to calculate the 
expected number of wide binaries as a function of angular sepa- 
ration. To this end, the above-described WW-technique requires 
the modification that we address now. 

Let ip(6)d9 be the number of wide binaries observed with an 
angular separation between 9 and 6 + d6. For a given distance D, 
we have ip(6)d6 — tfr(s)ds and s — DO. (Note that our unit for s 
and D is pc, whereas 6 is in rad.) The latter expression can be 
used to write the reduced separation distribution Q(s) as Q(D9). 
This introduces an explicit dependency on D, so we need to in- 
corporate Q(Dff) into the integration over D in the formula for 
the effective volume (Eq.l26T). This incorporation is the reason it 
is no longer possible to separate the reduced separation distribu- 
tion from the effective volume in a formal way like in Eq.[T"5l 

Furthermore, we need to modify the limits in the integra- 
tion over D in Eq. [15] Recalling that for a given semi-major 
axis a, the average observed projected separation is (s) = 0.98a 
(Eq. [22b . it appears to be appropria te to let the integ ration lim- 
its run from (s)min/# to (s) m&x /6 (lGarnavichlll99ll) . We take 
(s} min = 0.98a rain = 9.8 ■ 10~ 4 pc and <s) max = 0.98a T (for the 
calculation of the tidal limit aj see AppendixIBl. 

Putting it all together, we may write the number of observed 
wide binaries as a function of angular separation as 



« WB fl j dDD 3 p(D)Q(De)jjdM 1 dM2®(M l )®(M2)A43) 

(S)min/e M min (D) 



where we have included an additional factor D into the integra- 
tion over D, because ds = Dd9 for a given D. 

The model-2PCF is now determined by adding the number 
of physical pairs, calculated by Eq. |43] to the number of pairs 
expected from a random sample given by Eq. [6] 



m + Pi® 

P(9) 



1 



(44) 



The two free parameters - «wb and A - are determined in a least- 
square sense by fitting the measured 2PCF w con -(ff) to the model- 
2PCF w mo d(#) as described more in detail in the next section. 

4.4. Fitting procedure 

We determine the two free parameters of the model, «wb and 
A, by means of a Levenberg-Marquardt nonlinear least-square 



set of structure 




A 


Q 


x 2 /y 


/\VB 


parameters 


[pc- 3 ] 








% 


set 1 


0.0052 


1.00 


0.29 


1.14 


10.2 


set 2 


0.0069 


1.01 


0.30 


1.13 


13.6 


set 3 


0.0061 


1.01 


0.29 


1.14 


12.1 



algorithm (Lourakis 20040, which minimises the value of the 
X 2 defined as 



N 



: \ i / w con (6i) - w m0 d(6i) 



(45) 



The data is binned in steps of A8 = 6i+\ - 0, = 1", where 
Q\ - 6mm = 2" and On = #max = 30" as imposed by the reso- 
lution limit of the SDSS and maximal distance in the Neighbors 
table, respectively. As a result, we have N = 28 here. Given the 
approximative character of our study, the use of the this standard 
definition of the \ 2 is appropriate, even though the values of the 
2PCF at different angular separations are not strictly independent 

of each other. 

Following Press et al] dl992l) . we use the incomplete gamma 
function Q(x 2 \v) with v - N - 2 - 26 degrees of freedom as 
a quantitative measure of the goodness-of-fit. (N-2 because the 
model has two free parameters: «wb and A..) Values of Q near 
unity indicate that the model adequately represents the data. 

4.5. Confidence intervals 

To estimate the uncertainties of our best-fit valu es, we use 
Mont e Carlo confidence intervals (MCCRs) (e.g. iPress et all 
1992, §15.6). Assuming Poissonian errors, we generate 10000 
synthetic data sets by drawing the number of unique pairs "ob- 
served" in the Ai-th synthetic sample F from a Poisson dis- 
tribution with mean F(0), where F(6) is the observed number 
of pairs, not corrected for edge effects due to survey holes. The 
"true" number of pairs in a synthetic sample is then given by 
dividing by F™ (see 

For each synthetic sample, we determine best-fit values for 
the model's free parameters, nLL and A^ k \ in a least-square sense 
as described above. A /7%-MCCR is defined by the line of con- 
stant^- 2 which encloses p% of the best-fit values in the «wb ver- 
sus A plane (see Fig. [6j. The confidence intervals of «wb and A 
are then given by the orthogonal projection of the MCCR onto 
the corresponding axis. 

We include in our error estimate only the uncertainties stem- 
ming from the pair counts in F{9). Neither the uncertainties in 
the Galactic structure parameters nor those in the LF are taken 
into account. 



5. Results 

5. 1 . Analysis of the total sample 

The results of the analysis for the total sample are shown in 
Fig. [5] and Table [3] Figure [5] shows the 2PCF estimate w con (9) 
and the CDD y{9). The statistical uncertainties calculated ac- 
cording to Eq.[13]are shown as vertical lines. A strong clustering 
signal out to at least 9 = 10" is evident, whereas the CDD sug- 
gest that there are pairs in excess of a random distribution up to 



We use levmar-2.2. 
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Table 2. Parameters of the final sample and of the subsamples 



(Sub)sample 


n 


I 


b 


Nobs 






/Halo 1 " 




[deg 2 ] 


[deg] 


[deg] 




[pc] 


[pc] 


% 


total 


668.58 


175.6 


81.6 


669 843 


1555 


1.07 


30.0 


r < 20.0 mag 


668.58 


175.6 


81.6 


535 595 


1380 


1.06 


25.2 


r < 19.5 mag 


668.58 


175.6 


81.6 


425 674 


1235 


1.06 


20.6 


left 


334.57 


191.2 


73.7 


326 333 


1495 


1.08 


28.3 


right 


334.01 


102.5 


84.8 


343 510 


1625 


1.05 


31.0 


A 


78.94 


180.1 


68.3 


77 008 


1460 


1.09 


26.8 


B 


78.10 


167.5 


75.3 


76 011 


1510 


1.08 


28.6 


C 


78.99 


136.0 


79.8 


76 686 


1565 


1.06 


29.9 


D 


78.81 


93.3 


78.4 


82 780 


1635 


1.05 


30.6 


E 


88.23 


208.0 


69.7 


87 751 


1485 


1.08 


27.4 


F 


88.30 


210.6 


78.6 


85 563 


1545 


1.07 


29.5 


G 


87.95 


215.2 


87.5 


88 440 


1615 


1.05 


31.0 


H 


88.27 


33.4 


83.6 


95 604 


1700 


1.04 


31.8 



a Calculated using the Galactic structure parameter set 2. 
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Projected angular separation 9 in arcseconds 

Fig. 5. 2PCF as inferred from the total sample (solid circles, left 
ordinate) and the corresponding CDD (open circles, right ordi- 
nate). Poisson errors are indicated as vertical lines. Model curves 
for the three Galactic structure parameter sets are plotted, too, 
but as the differences between them are marginal they lie one 
upon the other, giving a single solid line. 



maximum angular separation examined, that is, up to 30". The 
outlier at 8 — 9" is probably a random fluctuation. We also plot 
in Fig. [5] best-fitting models using the three Galactic structure 
parameter sets described in Sect. 14.2.11 

The best-fit values of the two free parameters, «wb and A, are 
tabulated in Table [3] for the three Galactic structure parameter 
sets. The power-law index A appears to be quite independent of 
the set we choose. The wide binary density «wb, on the other 
hand, shows some variation. The difference between the sets 1 
and 2 reflects, for the most part, the difference in the overall 
normalisation p(0) of the density distribution, whereas in set 3 
the unequal halo normalisation with respect to the other two sets 
also contributes to the difference in «wb- 



1.2 




4.5 4.7 4.9 5.1 5.3 5.5 5.7 5.9 



Total wide binary number density n wB in (10 pc) 3 

Fig. 6. Distribution of the best-fit values for the total sample us- 
ing the structure parameter set 2 obtained by the Monte Carlo 
procedure described in the text. The solid contours (MCCRs) 
are lines of constant^ 2 and enclose 68.3%, 95.4%, and 99.7% 
of the best-fit values. 

We also list the goodness-of-fit Q and the corresponding re- 
duced chi-square values (x 1 divided by the degrees of freedom 
y) in Table [3] All three sets of Galactic structure parameters 
give equally good fits, whereas the reduced chi-square values, 
which are only slightly more than unity, indicate that we have 
not severely underestimated the uncertainty in the 2PCF. 

In Fig. [6] we show, representative of the other structure pa- 
rameter sets, the distribution of the best-fit values from the syn- 
thetic samples using set 2, our standard set. In the same figure 
the 68.3% (lcf), 95.4% (2a-), and 99.7% (3a) MCCRs are also 
plotted. Quoting 95.4% confidence intervals throughout, we find 
for our final sample using the Galactic structure parameter set 2 

«wb = 0.0052!^ pc~ 3 and A = 1.00!£g. (46) 

The power- law in dex A of semi-major axis distribution is con- 
sistent with lOpikl 's law (A = 1) up to the Galactic tidal limit, 
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Fig. 7. MCCRs for the three Galactic structure parameter sets 
corresponding to the total sample. Long dashed line: set 1; solid 
line: set 2; short dashed line: set 3. The crosses indicate the 
best-fit values inferred from the true observed number of pairs, 



whereas the number density «wb corresponds to a wide binary 
fraction with respect to all stars (no colour-cut) in the solar 
neighbourhood of 



_ 2»WB _ .„ 9 +1.2<w 

/wb = - 1 ( J.2_ 10 %, 



(47) 



where «* ^ 0.10 pc 3 is the t otal local stellar n umber density, 
i.e. the integral over the whole JahreiB & Wielen! LF. (The local 
wide binary density «wb corresponds only to wide binaries with 
both components having g - r > 0.5 mag.) 

We show the confidence regions corresponding to the three 
sets of Galactic structure parameters in Fig. [7] The confidence 
intervals for A agree for all the sets, which demonstrates that we 
can determine A and its confidence intervals reliably - provided 
that no systematic error has crept into our analysis. The wide 
binary density «wb, on the other hand, is more sensitive to the 
exact values of the Galactic structure parameters than A is, so the 
values we have derived for the wide binary density «wb and the 
fraction /wb should be viewed with caution. We are, however, 
confident that the true value of the wide binary density «wb is 
within a factor of 2 of our derived value. 

Since the structure parameter set 2 is - as far as known to the 
authors - the only one systematically corrected for unresolved 
multiplicity, we give that set m ore weight. From n ow on, we use 
the bias-corrected values from Jur ic et al . (2008), i.e. our stan- 
dard set 2, exclusively. 

Having determined the two free parameter in our model, we 
can calculate the number of wide binaries in our total sample as 
a function of the projected separation s. We show the distribu- 
tion in Fig. |8]whose shape is largely dominated by the effective 
volume, i.e. by selection effects. In particular, this distribution 
implies that we have observed 830^15 very wide binaries with a 
projected separation larger than 0. 1 pc in our total sample, while 
none are expected to be found beyond 0.8 pc, given our selection 
criteria. However, that extremely wide binaries with projected 
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Fig. 8. Distribution of wide binaries as a function of projected 
separation s that is expected to be observed in the total sample, 
given the model assumptions and the selection criteria we have 
adopted (see text). 
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Fig. 9. 2PCF and CDD for the r < 20.0 mag and r < 19.5 
mag subsamples together with the corresponding model curves. 
Up-pointing triangles (long-dashed lines) are for the r < 20.0 
mag subsample, whereas down-pointing triangles (short-dashed 
lines) are for the r < 19.5 mag subsample. The symbols were 
shifted apart by 0.25" for better visibility. 



separations of more than 1 pc can exist in the Galactic halo has 
recently been confirmed bv lOuinn et alj d2009l) . 

5.2. Differentiation in terms of apparent magnitude 

To test the consistency of the standard model, we used two sub- 
samples having brighter upper apparent magnitude limits than 
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Fig. 10. MCCRs for the r < 20.0 mag (long dashed contours) 
and r < 19.5 mag (short-dashed contours) subsamples. For com- 
parison, the MCCR for our final sample is also shown (solid- 
contours). The crosses have the same meaning as in Fig. [7] 



Table 4. Best-fit values: r < 20.0 mag and r < 19.5 mag 
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Fig. 11. 2PCF and CDD for the left and the right subsamples, 
together with the corresponding model curves. Left-pointing tri- 
angles (long-dashed lines) are for the left subsample, whereas 
right-pointing triangles (short-dashed lines) are for the right sub- 
sample. The symbols were shifted apart by 0.25" for better visi- 
bility. 



Subsample 



[pc~ 3 ] 



Q 



x h h 



r< 20.0 mag 0.0057 1.15 0.18 1.25 11.3 
r< 19.5 mag 0.0064 1.22 0.34 1.09 12.7 



Table 5. Best-fit values: left and right 



Subsample n WB 
[per 3 ] 



Q 



x h /\VB 

% 



left 
right 



0.0055 
0.0046 



1.00 0.55 
1.03 0.82 



0.94 
0.74 



10.8 
9.1 



our total sample and analysed them in the same way as the to- 
tal sample itself. We set the upper apparent magnitude limit to 
r = 20.0 mag and r = 19.5 mag - all other selection criteria (in- 
cluding the lower apparent magnitude limits) remain unchanged. 
The main parameters of these subsamples are listed in Table [2] 

In Fig. |9] we show the 2PCF and the CCD as in Fig. B 
together with the corresponding model curves for standard set 
2 (see caption). The best-fit values are listed in Table |4] The 
MCCRs of the r < 20.0 mag and r < 19.5 mag subsamples 
are shown in Fig.fTOl 

If the model were self-consistent, we would expect that the 
best-fit values agree with each other within their uncertainties. 
Figure [10] indicates that this is not the case: It appears that the 
best-fit values are systematically shifted to higher densities and 
larger power-law indices when using a brighter upper magnitude 
limit. As we discuss in Sect. 16.21 this inconsistency is most likely 
an artefact of the (oversimplifying) model assumptions and does 
not entirely undermine our results. 
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Fig. 12. MCCRs for the left (long-dashed contours) and right 
(short-dashed contours) subsamples. For comparison, also the 
MCCR for our final sample is shown (solid contours). The 
crosses have the same meaning as in Fig. [7] 



5.3. Differentiation in terms of direction 



In a previous study, ISaarinen" & Gilmord rfl989h found that the 
binaries appear to be highly clumped in the NGR However, it has 
not become entirely clear whether this patchiness of the wide bi- 
nary distribution in the sky is a real physical characteristic of the 
wide binary population or if it is due to statistical fluctuations. In 
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Table 6. Best-fit values: A to H 



Subsample 


«WB 


A 


Q 


T~i 

x'h 


/WB 




r — In 

[pc 3 ] 








% 


A 


0.0050 


1.27 


0.92 


0.64 


9.9 


B 


0.0054 


1.00 


0.96 


0.57 


10.7 


C 


0.0041 


1.02 


0.36 


1.07 


8.2 


D 


0.0046 


1.19 


0.29 


1.13 


9.1 


E 


0.0077 


0.73 


0.24 


1.18 


15.2 


F 


0.0045 


1.14 


0.37 


1.07 


8.9 


G 


0.0054 


0.83 


0.67 


0.86 


10.6 


H 


0.0042 


1.27 


0.64 


0.88 


8.4 



principle, we can check whether the wide binary density varies 
with position in the sky by dividing our sample into subareas. 

To begin with, we divide our sample in two halves by cutting 
it along the a = 185° meridian. In the following, we refer to the 
subarea with a < 185° as the "left" subsample and the other half 
as the "right" subsample. The subsamples' main parameters are 
listed in Table [2] 

In Fig. QT|we show the 2PCFs inferred from the left and the 
right sample, as well as the corresponding CDDs. It appears that 
there are a few more pairs in excess of random in the left half. In 
the same figure also the corresponding model curves are plotted. 
The best-fit values are listed in Table [5] 

In Fig. [12]we show the MCCRs of the left and the right sub- 
samples. The left subsample indeed shows a higher wide binary 
density than the right one. The difference is significant at the 3<x 
level. The power-law indices, on the other hand, do agree in both 
subsamples. Regarding the total sample, the right half differs sig- 
nificantly (at 3cr), whereas the left half is more consistent with 
it. 

The difference in the wide binary density between the left 
and the right halves is probably a real feature, as any inadequa- 
cies of our model (e.g. inaccuracies of the stellar density dis- 
tribution we use) should affect both halves in almost the same 
manner. To examine this apparent positional dependency of the 
wide binary density in more detail, we divide our sample fur- 
ther into eight subsamples, each covering 10° x 10°. They are 
labelled from A ("upper left") to H ("lower right") as suggested 
in Fig. Q~3] where the abscissa can thought of representing the 
right ascension from 165° ("left") to 205° ("right") and simi- 
larly the ordinate represents the declination from 22° ("bottom") 
to 42° ("top"). The subsample's main parameters are again sum- 
marised in Table|2] 

Figure [13] shows the 2PCFs estimate and the CDDs with the 
best-fit model curves. The corresponding best-fit parameters are 
listed in Table [6] Some subtle differences are apparent between 
different subsamples. While the CDD in A, B, C, and G are re- 
produced well by the model, the CDD in D, F, and H appears to 
be too flat. In those subsamples it seems that almost no physical 
pairs are present at angular separation s over 15", in agreement 
with the findings of lSesar et all d2008l) (see Sect.lOTi. 

The subsample E appears to be an outlier, because it contains 
almost twice as many pairs in excess of random as the seven 
other subsamples. The listed best-fit values confirm the peculiar 
character of subsample E having the highest wide binary density 
of all eight subsamples, together with the lowest power-law in- 
dex. As far as the authors can judge, there is no obvious special 
feature (e.g. open star cluster) in subsample E that could cause 
this anomaly. 

The 2<t MCCR are shown in Fig. [14] Except for the outlier 
E, all the other subsamples are quite consistent with each other. 
Also, no obvious trend, e.g. with Galactic latitude b, is apparent. 
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Fig. 14. 95.4% MCCRs for the subsamples A to H. The crosses 
have the same meaning as in Fig. [7] 

Table 7. Best-fit values: excluding subsample E 



Subsample 




A 


Q 


3 , 

x~/ v 


/WB 




[pc- 3 ] 








% 


total\E 


0.0049 


1.06 


0.40 


1.04 


9.6 


left\E 


0.0049 


1.12 


0.76 


0.79 


9.8 



To what extent is the subsample E responsible for the difference 
between the right and the left subsamples? 

To answer this question, we repeated the analysis of the to- 
tal and left (sub)sample excluding subsample E. The results are 
listed in Table [7] The wide binary density drops significantly 
to a value more consistent with the right subsample. Thus, we 
conclude that the difference between the right and the left sub- 
samples is largely caused by the outlier E. Apart from subsample 
E, the wide binary densities in different directions appear to be 
consistent with each other. The reason for the high density in 
subsample E remains unclear. However, a statistical fluctuation 
cannot be ruled out to a level better than 2cr. 

6. Discussion 

6.1. External (in)consistencies 

In this section we compare our results with those of previous 
studies. Our first point concerns the observed angular 2PCF. 
From the CDD of our total sample shown in Fig. [5] we have 
noted pairs in e xcess of random up to 6> max = 30". This is 
contradictory to ISesar et a! ] (I2008h who find that there are es- 
sentially no physically bound pairs with an angular separation 
9 > 15". (We note, however, that this is in accord with some 
of our subsamples, see Fig. Qj]) We considered several possible 
reasons for this apparent discrepanc y, such as differences in the 
selection criteria of the Ses ar et al.l sample with respect to our 
sample, underestimation of the total area of holes in our sample, 
or an overcorrection of edge effects due to those holes, but we 
convinced ourselves that none of them can account for this dis- 
agreement. A real physical difference in the wide binary popula- 
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Fig. 13. 2PCF (solid circles) and CDD (open circles) for the subsamples A to H, together with the corresponding model curves. The 
subfigures are arranged as the corresponding subareas would be seen on the sky. 
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Fig. 15. 2PCF and CDD as inferred from our total sample as in 
Fig. H] Th e model curve is fo r the broken power-law model from 
Lepine & Bon giornol d2007l) (see text). 



tion studied, e.g. caused by the fact that lSesar et all s sample is 
largely disc-dominated, while our final sample has a substantial 
hal o contribution of ~ 30%, can be excluded, as previous stud- 
ies dLatham et al.ll2002tlCh"aname & Gouldl2 004) found that the 
disc and halo wide binary populations are reasonably consistent 

in their statistical characteristics. 

Regarding the wide binary fraction, ISesar et all also find a 
much smaller wide binary fraction of below 1 % (decreasing with 
height above the Galactic plane), as compared to our/wB ~ 10%. 
Even if all pairs with a semi-major axis larger than 3 000 AU (be- 
yond the "break") were ignored, we woul d still find a wide bi- 
nary fraction of 4.6%. On the other hand, Lepine & Bon giornol 



(2007) give a binary fraction of at least 9.5% for separations 
larger than 1 000 AU, which is in rough agreement with our re- 
sults (setting a min = 1 000 AU we find / WB » 8.3%). 

Studying a much smal ler regio n co ntaining br i ghter stars as 
compared to our sample, WW87] and iGarnavichl (Il988l Il99lh 
found an un physically large wide binary density in the direction 
of the NGP. ISaarinen & Gilmord (1 19891) attribute this overden- 
sity to a large statistical fluctuation. The region where this over- 
density was found would be located in our subarea H. But as 
we probe fainter stars, no density enhancement is evident in that 
subarea. The slight overdensity w e found in subarea E is differ- 
ent from that previously noted bv IWW87I and IGarnavichl 

As to the separation distribution, previous studie s hav e 
found that the observational data are described by Opik's 
law (A _ 1 ) up to a certa in maximum separation ("break"): 
ILepine &" Bongiornol (120071) find this break to be around 3 500 
AU , beyond which a steeper slope, with A » 1.6, should ap- 
plv. lSesar et alJ (l2008h basically agree with Lepine & Bon giornol 
and find in addition that the maxi mum separation increa ses with 
height above the Galactic pla ne. iPoveda & Alien! d2004l) divided 
their wide binary catalogue (IPoveda et alJ fl994) into two sub- 
samples consisting of the oldest and youngest systems, respec- 
tively, and find that the "maximum major semiaxis for which 
Opik's distribution holds is much larger for the youngest bi- 
naries (a m - 7862 AU) than for the oldest (a m = 2409 AU)". 
|Ch aname & G ould (2004) find the data to be described by a sin- 
gle powerlaw with an index of approximately 1.6 for their disc 
and halo samples. However, they note a "puzzling" flattening 
of the distribution of the disc bin aries between Wand 25". As 
already suggested bv lSesar et all this flat range might be the do- 
main where lOpikf s law is valid. I n view of the substantial un- 
certainties inherent in studies by ILepine & Bongiornol and by 
|Ch aname & Gouldi the two studi es appear to be broadly con- 
sistent. Similar to our study. IGarnavichl d 1 99 ll) assumed that the 
semi-major axis distribution is described by a single powerlaw. 
He finds the power-law index to be 0.7 + 0.2 for his NGP sam- 
ple covering nearly 240 square degree. At intermediate Galactic 
latitudes he finds the slope to be steeper: A — 1.3 + 0.2. His data 
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seem to favour a (somewhat unrealistic^ 4 ]) "cutoff" around 0.1 
pc, especially for the NGP sample. 

The question of a break in the separation distribution 
has much been discuss e d in the context of DM constraints. 
Wasserman & Weinberg d!991l) . for example, show that the ob- 
servational data suggest a break but they do not require one sta- 
tistically. Looking at the CDD for our total sample in Fig. [5] we 
note that our best-fit model slightly overestimates the number 
of pairs in excess of random at larger angular separations, and 
a more curved model line would be preferred by the data. This 
could indeed be interpreted as a hint that the semi-major axis 
distribution is broken, because having a steeper power-law index 
from a certain semi-major axis on would result in a flatter CDD 
model curve. In principle, we could easily fit a broken power- 
law to our data as well. However, too many free parameters will 
only destabilise our modelling. Instead we check whether our 
data are also consis tent with the specific broke n power-law dis- 
tribution found by Lepine & Bon giornol d2007l) . To this end we 
use a broken reduced semi-major axis distribution of the form 

q(a) — q\(a) + qi{a) — c^a ®(a c — a) + CA 2 a~' h -®(a — a c ) (48) 

with the Heaviside function that introduces a break at semi- 
major axis a c . It is normalised, as in $4.1.11 according to 

r*cij r*®c f° a "r 

I q(a)da — I q\{a)Aa + I qi(a) da = 1 . (49) 

Using the values given bv lLepine & Bongiornol d2007l) . A\ = 1, 
A2 = 1 .6, and a c = 0.02 pc, we find for the normalisation con- 
stants 



cm = 



and 



In a c — In a„ 



(ttr/flc) 2 ~ 1 
l-A 2 



ca 2 - cx.a 



a 2 -i 



(50) 



(51) 



The distribution of the projected separations Q(s) then splits into 
a sum, too: 



Q(s) 



• s 



-M 



C Al (s) pc 1 + 



pc 



C A2 (s) pc" 



(52) 



where CaX s ) is defined as in Eq. [24] with the limits chosen ap- 
propriately when integrating over 77. 

We now keep the parameter given by Lepine & Bong iornol 
fixed so that our broken power-law model has only one free pa- 
rameter: the wide binary density «wb- In Fig. [15] we show the 
best-fit model curves corresponding to «wb = 0.0075 pc~ 3 . The 
broken power-law model also gives a decent fit to the 2PCF in- 
ferred from our total sample. (One should not be misled by large 
discrepancy between the model and the CDD showing the cumu- 
lative difference between the 2PCF and the corresponding model 
curve.) With a goodness-of-fit of 2 - 0.0047 and a reduced chi- 
square of x 2 / v - 1-85, the broken power-law model provides a 
significantly worse bestfit to the data than our single power-law 
model. Nevertheless, as the errors in the 2PCF are slightly un- 
derestimated, the model is not put into question until Q < 0.001 
dPress et al.lll992l their §15), and we conclude that the broken 
power-law distribution found by Lepi ne & Bongiornol (2007) 
can not be rejected with confidence either. 



14 " Physically, one doesn't expect to see a sharp cutoff." (Weinberg 
1990) 



6.2. Internal (in)consistencies 

FigureQJJshows that the best-fitting values for A and «wb are in- 
consistent, within the adopted random errors, when the limiting 
magnitude is varied. Both parameters are systematically shifted 
by roughly 20% when the limiting magnitude is changed by 1 
mag. Further inconsistencies were met when we tried to differ- 
entiate our total sample with respect to colour (not shown here). 
This very likely means that there are unaccounted for systematic 
errors in the modelling, resting on oversimplified assumptions. 

A basic assumption on the properties of wide binaries was 
that both components are drawn at random from the same stel- 
lar LF. There is, however, some evidence that this is not quite 
correct: Goul d et al. I (fl995l) noticed that "binaries are bluer and 
more distant than one would expect if t hey were formed of 
random combinations of field stars". Also Lepine & Bon giornol 
d2007l) find "that the luminosity function of the secondaries is 
significantly different from that of the single stars field popula- 
tion, showing a relative deficiency in low- luminosity (8 < My < 
14) objects". Similarly. ISesar etal] d2008l) report that blue stars 
(g - i ^ 2.0 mag, corresponding roughly to g - r 5 1.4 mag) 
that are a member of a wide binary, have more blue companions 
than expected from the LF. (For red stars, however, they find that 
the components are draw n randomly from the LF.) Moreover, it 
has long been known (e.g. Bahcall & Soneirall 19801 their Fig. 2) 
that stars of early spectral type have smaller scale heights than 
late type stars; therefore, the assumption that the stellar density 
distribution depends only on distance, and not on the luminosity 
of the stars, might be an oversimplification as well. 

On the other hand, the internal inconsistency discussed is 
overemphasised in Fig. [10] because we only include the statisti- 
cal errors stemming from the pair counts when determining the 
MCCRs. If we also included the errors in the Galactic structure 
parameters and those from the LF, the MCCRs would be consid- 
erably larger, and the inconsistency evident from Fig. [TOl would 
not be as severe as it seems. That a variation in the Galactic 
structure parameter has a significant impact at least on «wb, is 
apparent from Fig. [7] The systematic drift towards larger A and 
"wb, when adopting a brighter upper apparent magnitude limit 
(cf. Fig. [Tol l, still points to some inconsistencies in our model. 
We think, however, that these inconsistencies do not entirely 
undermine our major results, although the quoted uncertainties 
might be considerably larger. 



7. Summary and conclusions 

We have derived the angular 2PCF for nearly 670000 SDSS 
stars brighter than r = 20.5 mag and redder than g - r = 0.5 
mag in a region of approximately 670 square degrees around the 
NGP. Various corrections had to be made for quasar contami- 
nation, survey holes, and bright stars. There is an unambiguous 
correlation signal on small scales up to 30". We modelled this 
signal by a modifie d WW-technique , closely following a previ- 
ous study of lGarnavTc h1 (fT988l[T99l that was based on a much 
smaller sample of stars. The modelling involved a number of 
parametrised distribution functions: the spatial density and LF 
of stars, as well as the separation distribution of bi naries. The 
Galact ic model used is based on the recent study by Juric et all 
d2008h and the stellar LF derived bv iJahreifi & Wielenl (11997^ 
For the wide binary semi-major axis distribution we assumed a 
single powerlaw. Furthermore, essential assumptions are: (1) bi- 
nary stars follow the same density distribution as single stars, (2) 
both components of a binary are randomly drawn from the single 
star luminosiy function. These assumptions allowed a significant 
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simplification of the model. The Galactic structure parameters 
were fixed (we explore three different sets), while the local wide 
binary number density «wb and the power-law index A of the 
semi-major axis distribution were left free; i.e., they have been 
determined by a least-squares fitting algorithm. 

The best fit to the observed angular 2PCF of the total sample 
was obtained with A » 1.0, which corresponds to the canoni- 
cal iQgjkl law, and «wb ~ 0.005 pc~ 3 , meaning an overall lo- 
cal wide binary fraction of about 10.0% in the projected separa- 
tion range of 0.00 1 pc (200 AU) t o 1 pc (Galactic tidal limit) 
Previous studies dPoveda & Allenl l2004t ILepine & Bongiornol 
120071: ISesar et all 12008) have also found the data to be consis- 
tent with Opik's law, but only to a maximal separation that is 
considerably smaller than the Galactic tidal limit. Beyond that 
maximum separation, the distribution continues with a steeper 
decline. We have shown that o ur data are also consi stent with 
the broken powerlaw found by Lepine & Bon giornol although 
with a fit of lower quality (though involving only one free pa- 
rameter, namely, «wb)- Given this ambiguity, we are not able 
to put any strong constraints on the presence of a break in the 
wide binary separation distribution, which is regarded as one of 
the most interesting aspects of wide binarie s. As to the wide bi- 
nary f raction, we are in good accord with ILepine & Bongiornol 
d2007l) . whereas an apparent discrepancy with lSesar et al.l {2008) 
remains. 

A differentiation of the sample with respect to limiting ap- 
parent magnitude turned up a systematic dependence of the bi- 
nary parameters on the sample depth, which is most probably 
an artefact caused by oversimplified model assumptions. This 
conjecture is strengthened by the impossibility obtaining self- 
consistent results when differentiating with respect to colour. We 
conclude that one or more of the simplifying assumptions put 
into the model (e.g. that the luminosities of the binary compo- 
nents are independent of each other and draw randomly from the 
single-star LF or that binaries' density distribution follows ex- 
actly the single-star density distribution and that the density dis- 
tribution is independent of the stars' luminosities) are not quite 
correct. 

Differentiating the sample in terms of direction in the sky did 
yield some modest but non-significant variations. Only in one 
direction (subarea E) was an unexplained overdensity found. 

While we have shown here that the stellar angular 2PCF, as a 
complement to common proper motion studies, basically works 
and remains a viable tool for the study of wide binaries, it has 
become clear that this method is severely limited by the need for 
- even more - complex modelling. To relax those simple model 
assumptions would unduly complicate the analysis much further 
and probably no longer yield unique solutions. Any more effi- 
cient progress will indeed have to involve distance information 
to discriminate against unwanted chance projections. We plan to 
include distance information in our future work. In spite of the 
limitations of the present, simple modelling of the angular 2PCF, 
we think that the general result for the total sample derived here, 
i.e. A x 1 (Opik law) and /wb ~ 10% among stars having a spec- 
tral type later than G5, to within an uncertainty of 10%-20%, is 
a robust result. 



Appendix A: Edge correction for holes 

To quantify the edge effect from holes in our sample, we assume 
that the holes are circular and flat (Euclidian) and that effects due 




Fig. A.l. Illustration of the geometry used to correct for edge ef- 
fects due to holes. The large circle represents the bounding circle 
of a hole or a bright star mask. A star lying close to the hole is 
indicated. The region within a distance 9 of that star (highlighted 
in grey) is partially inside the hole (dark grey intersection). 

to intersecting holes are negligible. Let A k (9) be the solid angle 
of the annulus of width 20 around a hole of radius k 

A k (0) = 4n9(e k + e) , with £ = {SH,BS), (A.l) 

where k stands for hole masks (survey holes: SH) or bright star 
masks (BS). Due to the holes in our sample we observe only a 
fraction F l ^(0) of all stellar pairs in the residual solid angle £2 
sparated by an angular distance 9 

F«(0)=1-^- 2 N k A k (0)(l-F H (O;9 k )), (A.2) 

<:=(SH,BS) 

where Fn(9; 9 k ) is the fraction of all pairs separated by 9 we 
observe in A k {9) 

2n r h+e 

F H (9;9 k )= = xf(x;9,9 k )dx. (A3) 

A k (9) Je k 

Here, / is the fraction of the area n9 2 of the disc with radius 9 
around a star that lies outside the hole (Fig. lA. lb 



fix; e,0k) 



n out (*;0A) 



1 



fiin(*;0,0*) 



(A.4) 



tt9 2 n6 2 

whereas Q; n = nO 2 - Q out is the area of that disc within the hole, 
i.e. the intersection. With the nomenclature given in Fig. lA. II we 
find 

Q oat = e 2 k a-yp + e 2 /3-yq, (A.5) 
where 



a = arcsin ( — I ; y8 = arcsin(-| 
[9,]' p \e) 



and 



p= p k -y 2 , q= ^ 2 -y 2 - 

Solving x = p + q for y using Maple gives 

— Jl9 2 h x 2 - x 4 
2x > * 



y 



+ 29 2 x 2 + 29 2 9 2 - e* . 

k k 



(A.6) 



(A.7) 



(A.8) 



where the integration in (IA.3I ) was performed with Maple, too. 
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Appendix B: Galactic tidal limit 

A crude estimate of the Galactic tidal limit aj, i.e. of the max- 
imum semi-major axis a max of a binary star with total mass M 
orbiting in the Galactic tidal field is provided by the Jacobi limit 
rj (sometim es also referred to as zero-v elocity surface or Roche 
sphere) (e.g. lBinnev & Tre maine 2008, §8.3) 



M 



n 



x3M(<D) 



^ D, 



(B.l) 



where D is the galactocentric distance of the binary system and 
Al(< D) is the Galactic mass enclosed within radius D. Then, 
Al(< D) is given by 



M(<D) 



-D 1.1 ■ 10' 



(IT 

P c , 



Mo, 



(B.2) 



where we have adopted the canonical value of the circular speed 
v c = 220 kms -1 . 

For a rough estimate of this distance range, it is useful to de- 
fine the effective depth of the sample as the most likely dis tance 
of an arbitrarly chosen star (Weinberg & Wasserman 1988). In 
a magnitude-limited sample the total number of stars expected 
to be observed in a (heliocentric) distance range between D and 
D + dD is 



M m . a (D) 



dN(D) = Q.dDD 2 p(D) 



I 



dM<D(M) . 



(B.3) 



M mm (D) 



where Q denotes the solid angle covered by our sample, p(D) 
the density distribution, and O(M) the stellar LF. The absolute 
magnitudes M mm (D) and M max (D) are given by Eq. 1771 and [28l 
respectively. The effective depth D P ff is then given by the media n 
value of the distribution AN{D) (Weinberg & Wasserman 1988) 



J^ s dN(D) _ fofdN(D) 
f °°dN(D) ~ g-dN(D) 



(B.4) 



where D m i n and Z) max are used for the numerical integration and 
are chosen to bracket the theoretical distance range probed by 
our sample, which - neglecting interstellar dust extinction - 
is readily found from the apparent magnitude limits, m^ m and 
wimax, an d me absolute magnitudes limits, M m \ a and M max , given 
by <3>(M). We take them to be 



D • = 10 ( " 

^min — rw 



-M m 



;+5)/5 _ 10 (15-21.5+5)/5 _ q g 



pc 



(B.5) 



and 



10 ( " 



,5)/5 = 10 (20.5-(-0.5) + 5)/5 a m kpQ ( g £) 



In Table|2]we list the values of D e ff for our various (sub)samples 
calculated using the Galactic structure parameter set 2. With the 
notation of 34.2.11 the effective galactocentric distance D e g can 
be expressed as 



D eS = r + D eS cos b — 2ro-D e ff cos b cos ( 
+zl + 2z £>eff sm b + D 1 ^ sin 2 b . 



(B.7) 



Its values are typically somewhat above 8 kpc. 

The last parameter we need to specify in order to estimate 
the Galactic tidal limit aj, is the total mass of the binary star 



M. Given the magnitude limits of our sample we may write the 
average mass of a star in the sample as 



j m£(m)dm J m(M)<D(M)dM 



<m> = 



(B.8) 



J £(m)dm J (D(M)dM 



where £(m) denotes the mass function and m(M) th e mass - 
luminosity relation which we take from iKroupa et al.1 dl993l) . 
The transformation into r magnitudes was performed as in 

£221 

Statistically, we expect the total mass of a binary star to be 
M = 2(m), and we finally estimate the Galactic tidal limit ac- 
cording to 



fl T (£> eff ,M) =! 



M 



3M(<D eff ) 



Oeff 



3.1 • 10 



-3 



iuy 


L>e.ft 


\Me) 


{ P c J 



pc. 



(B.9) 



The values of aj are around 1 pc and we list them for various 
(sub)samples in Table|2] 
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